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An investigation of the mesoscopic dynamics of chemical systems whose mass 
action equation gives rise to a deterministic chaotic attractor is carried out. A 
reactive lattice-gas model for the three-variable autocatalator is used to provide a 
mesoscopic description of the dynamics. The global and local dynamics is studied 
as a function of system size and diffusion coefficient. When the diffusion length is 
I comparable to the system size phase coherence is maintained but the amplitudes 

of the oscillations are uncorrelated due to interaction between fluctuations and the 
instability of the chaotic dynamics. If the diffusion length is small compared to the 
system size then phase turbulence serves to destroy the noisy global attractor. 
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I. INTRODUCTION 
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■ The dynamics one observes on macroscopic scales has its origin in the molecular collision events 

CN ■ that occur on the microscopic level. Under far-from-equilibrium conditions even the macroscopic 

dynamics can take complex forms and show chaotic or turbulent behavior. Macroscopic chaotic or 
ordered structures arise because the correlations extend to macroscopic scales. A complete under- 
standing or description of the dynamics entails an examination of the system on the microscopic, 
C3 ' or at least mesoscopic, length scales where internal fluctuations enter naturally in the dynamical 

description. 

When the macroscopic deterministic dynamics is chaotic it is perhaps even more intriguing to 
O ' consider how such structured chaotic motion for the macroscopic fields arises from the apparently 

, random collision processes at smaller scales. We are then led to consider the interactions of 

fluctuations on small scales with the intrinsic instability of the chaotic dynamics on the macroscopic 
scales. 

In this article we use a reactive lattice-gas automaton to provide a mesoscopic description of 
k>( \ the reaction dynamics. We investigate a specific mass action chemical scheme, the three-variable 

j_j ■ autocatalator, which shows a period-doubling cascade to a chaotic attractor. The reaction scheme 

I and its bifurcation structure are described in Sec. |l|. The lattice-gas model naturally incorporates 

local particle number fluctuations as a result of random local reactive events and particle motion. 



Section [II provides details of the construction of the automaton dynamics and demonstrates that 



its mean field limit is the mass action rate law. The spatiotemporal dynamics in the chaotic regime 



is the subject of Sec. IV. The roles of internal fluctuations, diffusion and system size in determining 
the global and local dynamics are explored in this section. Section ^ presents a discussion of the 



automaton space and time scales and their relation to those in physical systems, while Sec. VI 
contains the conclusions of the study. 

The automaton dynamics considered here not only provides a more fundamental description 
of far-from-equilibrium, chaotic, chemical systems but also constitutes a scheme for the study of 
such systems in the mesoscale domain where standard reaction-diffusion equation models may be 
inappropriate. 



II. THREE- VARIABLE AUTOCATALATOR DYNAMICS 

In order to investigate the mesoscopic dynamics of systems exhibiting deterministic chaos it is 
necessary to consider reaction mechanisms with at least three species, apart from pool chemicals 
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whose concentrations are maintained at fixed values by external feeds and constrain the system 
to lie far from equilibrium. It is also important that chemical concentrations remain positive 
and follow mass-action kinetics]^ In this article we consider the three-variable autocatalator |2j 
which has a rich phase space structure including periodic and chaotic attractors^ and allows us 
to investigate certain aspects of the dynamics without unnecessary complications arising from 
bistability in the relevant parameter range.^ 

The three- variable autocatalator reaction mechanism is as follows 01: 
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The chemical dynamics is described by the variations in the U , V and W concentrations, while 
concentrations of A and B are taken to be fixed by feeds of these species; consequently thier 
concentrations can be regarded as bifurcation parameters. The mass action rate law corresponding 
to this mechanism is 

K2PuPI + I^ZPu - K4Pv , (2) 
I^APv ~ l^bPw ■ 

Here Pr(i) is the concentration of species t = u,v,w. The Ki, i = l,---,5, are effective rate 
constants that contain the concentrations of pool species. The system (||) has a steady state 

{pu,Pv,Pw) = (Ws,Ws, Ws) where Ms = liiVs/{K2vl + K3), Vs = K.QK^/Ki{K5 - Kl) and = K4Vs/k5. 

This steady state may be destabilized by a Hopf bifurcation at which a small-amplitude, oscillatory 
state emerges. In this study we choose K2 as the bifurcation parameter while the values of the 
other parameters are fixed at kq ~ 1.666667, ki = 0.153, K3 = 0.02, K4 — 4.0 and K5 — 1.0. With 
this choice of system parameters the Hopf bifurcation point occurs at K2 — — 15.68. The steady 
state is stable for K2 > k^- When K2 is decreased a period doubling cascade leading to chaos is 
observed (cf. Fig. |l|). 
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III. REACTIVE LATTICE-GAS MODEL 



The mesoscale description of the autocatalator is based on a reactive lattice-gas model of the 
dynamics. The method of construction of such models along with specific applications have been 
described earlier [^,|| and, consequently, only sketch of those features of the model that pertain to 
the present application will be given. 



^One of the first models of this type to be constructed is that due to Willamowski and Rossler, Ref. [Q 
^This model is an extension of the two-variable autocatalator developed by Gray and Scott, Ref. |^ 
^In the Willamowski-Rossler (WR) model the chaotic attractor coexists with a stable fixed point which 
can lead to interesting noise-induced transition processes but whose description can complicate the inter- 
pretaion of internal noise efi'ects on chaotic dynamics. Studies of noise-induced transition phenomena in 
the WR model can be found in Ref. Hi 
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Imagine that each chemical species Xr [U — Xi, V = X2, W — X3) resides on its own lattice 
£7. . Each node of the species lattice Ct contains channels corresponding to unit velocities directed 
along the links of the lattice from a node to its nearest neighbors.^ The nodes of the species lattices 
have identical node labels r so that several different species may reside at the same spatial point in 
the system. There is no exclusion principle in the automaton model used in this study; i.e., each 
velocity direction of a node on Cr may be occupied by any number of molecules qit- > 0.0 

Particle diffusion occurs through successive applications of deterministic propagation of particles 
in directions determined by their velocities to neighboring lattice nodes, and velocity randomization 
where velocities are shuffled and particles are assigned new velocities at random. We denote the 
combination of propagation and velocity reandomization on Ct by Vr, the diffusion operator. 

Consider a lattice C with Af nodes labeled (ri,r2, • ■ • ,rjv) and suppose many applications of 
the propagation and velocity randomization steps V have been carried out. When there is no 
restriction on the number of particles at a node on a lattice, the equilibrium distribution of n 
particles on Af nodes is an occupancy problem. The probability of obtaining occupancy numbers 
Q!(ri), a(r2), • • • , a{r^f) (a(ri) + a(r2) H h a{rj^) = n) on TV nodes is given by |^,^ 

n' 

P(a(ri),a(r2),---,a(rA^)) = ' -.A/"-". (3) 

Q;(ri)!a(r2)! • --a^Vj^y. 

The multimomial distribution (^) can be generated by the multinomial theorem 

(Pi+P2H VPN-T= X! -P(a(ri),a(i"2), • • • ,a(rA/-)) , (4) 

a(ri),Q(r2),---,Q(rAr) 

with pi — P2 = ■ ■ ■ — Pj\f — 1 /Af- To obtain the particle number distribution on a single node we 
sum over all the probabilities for the remaining Af — 1 occupancy numbers: 

P(a(ri))= P(«(ri),a(r2),---,a(r^)). (5) 

Q(r2),---,a(rAr) 

The constraint on occupancy numbers a(r2), ■ ■ • , a(rjv') in the sum is a{r2) + - ■ ■+a(rj\f) = n—a(ri) 
with < a{ri) < a{rj^). Therefore we obtain 

^ \ a(ri) / ^ n-Q(ri) 



where we have used the fact that 

V {n-a(ri))! ^„(„_(r,)) = [(_^ _ l)/A/]"-"(rO. (7) 

a{T2),■■■,a(r^r) 

Since diffusion is carried out independently on all species lattices Ct, we may simply append the 
label T to the node occupancy numbers and suppress the node label since all nodes are equivalent. 
Thus the particle number distribution at a node on Cr is binomial: 



*One may also include zero velocity or stop particle channels in the node description (cf. Ref. g). 
bee, Refs. i and §. In the actual implementation of the model the occupancy of a node is restricted 
to some large number of particles, much larger than the average occupancy of a node. This insures that 
the elements of the reaction probability matrix always lie between zero and one. 
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where n,- is the total number of particles on Lr- If we let A/" — > cxd and — > od in such a way 
that the ratio Pj = tLt/M remains finite, the binomial distribution P\,{ar) tends to the Poisson 
distribution 

i^p(«r) = ^^exp(-p.). (9) 

The automaton reactive dynamics at a node is controlled by the reaction operator C whose ac- 
tion is encoded in the probability matrix P{a.\l3). Each element of P{a.\(3) gives the transition 
probability from a reactant particle configuration a = (ai, a2, ■ • ■ , ckj) to a product particle con- 
figuration (3 = (/3i,/32, • • • ,/?s), where s is the number of species. If reaction is a rare event and 
diffusion is able to maintain spatial homogeneity over the lattice the local particle distribution will 
be homogeneous (binomial or Poisson). In this case the (discrete) mean field rate equation for the 
concentration changes is easily written in terms of the reaction probability matrix as 

p,(^ + l)-p,(0 -^^(/3.-a.)^(a|m(a,p(i)) , (10) 

where P\^{oi, p{t)^ = ^TPhi'^T)\p=p^(t) with Pi^la.^) — Pb(^) or Ppictr)- In the limit n,- — > oo (Vr) 
and small reaction probabilities the mean field equation (|lO|) should agree with the phenomenolog- 
ical mass action rate law (||). The reaction probability matrix P(a|/3) can be constructed so that 
this mean field limit is obtained. If we restrict particle number changes at a node to increases or 
decreases by one particle on each species lattice in accord with the mechanism of the reaction, one 
of the possible choices for the reaction probability matrix is 

P{ct\l3) = pi(Q:)5/3„^Q„-|-l<5/3„,a„<5/3„,a„ +P2(Q:)'5/3„,q„-1<5/3„,q„ + i5/3„,q„ 
+ P3(a)^/3„,a„<5/3„,Q„-l<5/3„,a„-|-l +P4(a)^/3„,a„i5/3„,a„5/3„,a„-l , 

where 

Pi{a) = h{Ko + Kia^) , 

P2(a) h[K2auayiay - 1) + K3a«] , 
P3{a) = hmay , 

P4,{ol) = hn^aw . (12) 

Here /i is a factor that determines the time scale in the automaton simulation. If (|l^) is inserted 
into (|l^) the mass action rate law (|^) is recovered. Thus we may now investigate the full automaton 
dynamics and compare its behavior to the autocatalator mass action rate law and reaction-diffusion 
equation.^ 

In summary, the automaton dynamics is specified by the compositions of the diffusion Dr and 
reaction C opertors: 

T = ' ' 



n(^r)^°C"=. (13) 



In this equation l-r and Ic are non-negative integers, indicating that the operations are repeated lr 
and Ic times, respectively. The reactive lattice gas model is a synchronously-updated, probabilistic 
cellular automaton. 



The reactive lattice-gas automaton bears some relation to birth-death master equation methods. See, 
for instance, Ref. jlo| 
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The simulations reported in this paper were carried out on triangular lattices with M = N x 
N nodes with periodic boundary conditions. On such triangular lattices the application of the 
diffusion operator V yields a diffusion coefficient D — 1/4 in units of lattice spacing squared per 
time step. For the composition rule in the diffusion coefficient of species r is given by = 
{lr/lc)D. In the simulations we have taken the diffusion coefficients to be equal: = Dy — D^j. 
The time scale factor was taken to he h = 0.001 . 



IV. SPATIOTEMPORAL DYNAMICS IN THE CHAOTIC REGIME 



A complete understanding of the the dynamics of chaotic systems on mesoscopic scales involves 
the consideration of a number of different aspects. For systems which are sufficiently large and 
contain many particles but are well-mixed, either through diffusion or mechanical means, one 
expects that the chaotic attractor will resemble that of the mass action rate law, apart from 
modifications of the attractor structure on small phase-space scales, [^t]-^ If the system is large 
and molecular diffusion is the sole mechanism responsible for eliminating concentration gradients in 
the system then the full spatiotemporal dynamics must be considered. Even in such a circumstance 
there is interesting information in the global (spatially averaged) concentration fields and the 
attractors corresponding to these global variables; for example, for coupled map lattices composed 
of chaotic maps non-trivial, low-dimensional, global attractors have been observed. The nature 
of such global attractors has also been studied in lattices of diffusively coupled Rossler oscillators 
where the local Rossler ODE dynamics is chaotic. The description of all aspects of the noisy 
mesoscale dynamics involves a study of the spatial structure that underlies the temporal evolution. 
This section is devoted to an investigation of some of these issues for the autocatalator where we 
focus on the dynamics for K2 — 11.52, a parameter value lying in the chaotic regime. 



A. Noisy chaotic attractor: global dynamics 

We observed above that if the system is small enough so that diffusion is able to maintain spatial 
homogeneity over the entire lattice but contains a sufficient number of particles that fluctuations 
in the mean concentrations are small, the automaton attractor constructed from the globally- 
averaged concentrations will closely resemble that determined from a solution of the mass-action 
ODE. Detailed studies of this type have been carried out for the Willamowski-Rossler model as a 
function of system size where the automaton rule has been modified to guarantee perfect mixing. 
[ p2t These calculations have confirmed that the gross structure of the chaotic attractor survives 
under the infiuence of internal noise, although noise can lead to destruction of fine-scale structure 
(e.g., band merging in banded attractors). There has been considerable discussion in the literature 
about the effects of internal noise in chaotic systems and detailed discussions can be found in Refs. 



As an example, for a lattice of size 32 x 32 nodes we show that the noisy attractor (cf. Fig. ^(b)) 
constructed from the spatially-averaged concentrations, 

p,(t) =iV-^ ^ p,(r,t) , (14) 

resembles the deterministic attractor shown in panel (a) of Fig. H, although there are noticeable 
effects due to fluctuations; as expected, the fluctuations tend to spread the dynamics in the unstable 
directions and obliterate some of the attractor structure. These effects diminsish as the system 
size increases provided spatial homogeneity can be maintained. Jl^ 

Next, we consider the attractor structure determined from the spatially-averaged concentrations 
as a function of the system size for a fixed value of the diffusion coefficient, Dr — 1/4. The results 
of a series of simulations on a set of x triangular lattices, with N = 64, 128, 256 and 512, are 
shown in Fig. 0. One observes that as N increases the global attractor increasingly resembles a 
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noisy periodic attractor whose width shrinks as the system size increases, at least on the time scale 
of the simulations. In the remainder of this section we examine the structure of the spatiotemporal 
dynamics that underlies this behavior. 



B. Spatiotemporal structure 

We shall now examine in some detail the spatiotemporal dynamics on a 128 x 128 lattice. The 
spatial structure that underlies global attractor discussed above is presented in Fig. ^. This figure 
shows the local concentration per node of species U for different times (indicated by open squares) 
along a trajectory segment of the global attractor given in Fig. |[ Strong concentration gradients 
develop during the fast parts of the cycle while only small gradients exist during the slow parts of 
the cycle. Of course, this behavior is anticipated in view of the ability of diffusion to homogenize 
the concentration variations that arise from local reaction. The slow part of the cycle, tg, comprises 
most of the average cycle time on the attractor, tc- For the parameter value under consideration, 
K2 — 11.52, we have tc ^ 2.3, while ts ~ 2.0 and thus for the fast part of the cycle tf ^ 0.3. 
Consequently, the diffusion lengths corresponding to these times are: ~ 32 and ^ 12 lattice 
sites, respectively. Since comparable concentration variations occur in the slow and fast parts of 
the cycle diffusion will clearly be able to more effectively homogenize the system in the slow part 
of the cycle than in the fast part of the cycle. 

The spatiotemporal structure for the same parameters but on a 512 x 512 lattice is shown in 
Fig. ||. Once again the strongest spatial inhomogeneities occur during the the fast part of the 
cycle but even in the slow parts, as is evident from the first three panels in the figure, one sees 
the small- amplitude inhomogeneities on the same length scale as in Fig. ^ more clearly since the 
system size is so much larger than the diffusion length. 

The spatial structure may be examined in more detail by considering the coarse-grained dynamics 
of the system. For this purpose we divide the system into cells with linear dimension L and define 
the coarse-grained concentration fields as, 

p,(r,t)=i-2 J2 pr{v,t), (15) 

re£x(f) 

where Ct{v) is a domain of the lattice centered on f containing L x L nodes. The time variations 
of Pu(f , t) for four cells with L = 32 are shown in Fig. |^. The chaotic oscillations in the individual 
cells remain in phase (at least over the time scale of the simulation]^ but their amplitudes vary 
erratically, especially during the transition from the slow part of the cycle to the fast part where 
fluctuations are strongly amplified. Such regions of the trajectory are expected to be most affected 

by noise arising from internal fluctuations .n 

Effects of this type were discussed earlier ||l9|] in a consideration of the global dynamics of coupled 
map lattices with chaotic or periodic local elements. In the chaotic regime any internal fluctuation 
(5/9T-(r,0) will grow roughly as 

<5p,(r,n,) = e^*^"^5p,(r,0) , (16) 

where ric — t/tc is the number of cycles on the chaotic attractor and A is the Lyapunov exponent. 
Here we have taken Pr (? , t) to be the coarse-grained density. If the system is well-stirred the particle 
fluctuations at each node of the lattice should be Poisson distributed as discussed in the previous 
section. For flnite diffusion, given the estimate of the diffusion length, one may expect Poisson- 
distributed fluctuations within the coarse-grained cells of size i = 32. The Poisson fluctuations 



^The total simulation time was 200 time units comprising about 100 cycle times on the attractor. 
^Discussions of the effects of external noise on different parts of the cysle for both chaotic and periodic 
dynamics can be found in Ref. [jL8| 
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within such coarse-grained regions provide a source of perturbations of the dynamics of magnitude 
Spr{f,0) = {prf'^/L which are then amphfied by the chaotic dynamics. One may then estimate 
the number of cycles on the attractor for the system to lose memory of its initial state as 

1 &p%{v)L 

where 5p*{r) is a value of Spriv^t) comparable to the "size" of the attractor in phase space. Using 
the cycle time on the attractor, tc = 2.3, and the fact that the Lyapunov exponemt is A = .229 
for K2 — 11.52 one finds that n* is of the order of a few cycles. Thus, given this time, one expects 
that regions separated by distances of order £c {2Dtcn*y^^ will have the amplitudes of their 
oscillations uncorrelated. One finds £c ~ 33 and thus the amplitudes in the 32 x 32 coarse-grained 
cells should vary randomly as observed in the simulations. If one then spatially averages over 
the coarse-grained cells one should see noisy periodic dynamics since the random amplitudes will 
average to a nearly constant value. Note that phase coherence has not yet been destroyed and the 
above arguments apply to the structure transverse to the attractor. 

It is interesting to compare the above dynamics with that in the periodic regime. The results 
for K2 = 12.8 are shown in Fig. |^ where panel (a) is the deterministic attractor and panel (b) is 
the attractor constructed from the global dynamics determined from a spatial average over the 
128 X 128 lattice. As expected, a noisy periodic attractor is obtained with strong amplification of 
the fluctuations in the fast part of the cycle. If one examines the coarse-grained dynamics that 
underlies the noisy periodic orbit (cf. Fig. ^ one sees local dynamics that is qualitatively similar 
to that for the chaotic dynamics described earlier: the internal fluctuations most strongly influence 
the the portion of the cycle where the dynamics changes from slow to fast, and phase coherence 
is maintained for this system size and value of the diffusion coefficient. In the periodic regime, 
from a linearized analysis of the dynamics, one expects the fluctuations about the deterministic 



motion to grow linearly with time. 11| The fact that the dynamics observed in our simulations 
is similar in both the chaotic and periodic regimes likey reflects the small magnitude of the positive 
Lyapunov exponent for the autocatalator chaotic attractor. This implies that rather large spatial 
regions may remain correlated even in the chaotic regime. Similar persistence of correlations due 
to small positive values of the Lyapunov exponent were noted by Brunnet et al. p5| in studies of 
coupled Rossler ODEs. 

In summary, the picture of the dynamics that is seen is consistent with the qualitative bifurcation 
diagram in Bohr et al. JTot ; for the finite-amplitude noise arising from internal fluctuations in the 
dynamics one sees a noisy periodic orbit in both the chaotic and periodic regimes. For the relatively 
small lattice size, 128 x 128, the system maintains phase coherence but noise interacts with the 
instability in the dynamics in directions along the unstable manifolds of the chaotic attractor to 
lead to erratic amplitude oscillations whose spatial average gives rise to a noisy periodic attractor. 



C. Destruction of phase coherence 

As the system size increases, or alternatively the diffusion coefficient decreases, it will become 
impossible for the system to maintain phase coherence over the entire spatial domain. [^,0 To 
investigate the passage to this regime we again consider 128 x 128 lattices but vary the ratio £r/£r 
to change the magnitude of the diffusion coefficient (and reactive time scale). In the automaton 
simulations we increase the value of £r keeping £r fixed so that the effective reaction rates are 
increased. 

The global attractors for = D/6 and D/10 are shown in Fig. |l^. They shrink to a noisy 
fixed point attractor as synchronization of the phase in space is lost. This is easily confirmed from 
an examination of the coarse-grained dynamics in Fig. now ones sees that both the amplitude 
and phase vary in the coarse-grained cells. 

If the description of the dynamics were based on the reaction-diffusion equation it would be 
possible to account for changes in the overall reaction rates and diffusion coefficients by scaling 
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space and time. If such scaling were valid, one could relate the automation simulations for fixed 
system parameters and varying system size to those for a fixed system size and varying parameters. 
This scaling is not valid as is evident from the results presented above. This signals a breakdown of 
the reaction-diffusion equation model for the fast reaction rates considered in this subsection. Here 
several reactions were carried out for each diffusion step giving rise to strong local fluctuations due 
to reactions which could not be smoothed by diffusion. 

If one considers dynamics in the periodic regime (k2 — 12.8 described above) and decreases the 
diffusion cocfhcient or increases the system size, destruction of phase coherence is observed. Such 
destruction of phase coherence on large scales is well known in oscillatory media and gives rise 
to phase turbulence. The present results for the chaotic regime show that even though the 
dynamics in small! correlated regions is locally chaotic, phase coherence on the chaotic attractor 
is not maintained on long scales, a result that parallels that in oscillatory media. Thus, our 
stochastic reactive lattice-gas simulations of the autocatalator do not show a robust global periodic 
attractor that survives as the system size increases. In certain parameter regimes, such robust 
global attractors were seen in deterministic systems of coupled Rossler ODEs. [|l^ 



V. SPACE AND TIME SCALES 



A number of aspects of fluctuations and their effects on chaotic dynamics which were described 
above depend on the space and time scales involved in the automaton simulations. In order to use 
these results to interpret the behavior of physical systems one must establish a link between the 
scales in the automaton simulations and those that arise in real systems. Any connection of this 
type must take into account the characteristics of a specific system so we shall simply present a 
few order-of-magnitude estimates for representative cases to illustrate the points. 

Consider a liquid-state (aqueous) reacting system where the reactants and products are small 
molecules. For a O.OIM solution, a volume of solution with a linear dimension of 0.01/i contains 
about three reactant molecules, a typical occupancy of an automaton node in our simulations. 
Thus, we may associate a linear dimension of £a — 0.01/i with the automaton lattice spacing. 
Taking the molecular diffusion coefficient to be _D ^ 10~^ cm^/s, the average time for a molecule 
to move by diffusion to a neighboring cell is £l/D = Ta ^ 0.1 /is. This time can be associated with 
an automaton time step. 

Next, we must consider the time scale of the reactive events in comparison with that for the 
non-reactive events, those that are largely responsible for the diffusion process, in order to see 
under what circumstances these scales may be applied. For this puropose we may use the average 
cycle time observed in the global dynamics. In the simulations presented in this paper the average 
cycle time was of order tc ~ 10'^ time steps, and this corresponds to tc ~ 0.1 ms, a rather fast 
chemical relaxation time for most typical oscillatory chemical reactions, but not unrealistic for 
reactions under certain conditions. 

More usually the cycle time of many laboratory examples of oscillatory reactions lie in the range 
1 — 100 s. One may then work backward and associate a time of 10~^ s with an automaton time 
step for a cycle time of tc = 1 s from which it follows that the automaton cell dimension is roughly 
10~^ cm if the same value of the diffusion coefhcient is used. In a cell volume of this size there 
are about lO'' reactant molecules so that the fluctuations in the automaton are about 10'^ times 
larger than in a system of this type. Of course, these estimates serve to confirm that in usual 
circumtances, on the space and time scales of many macroscopic experiments, internal fluctuations 
play only a minor role. However, they are crucial in homogeneous nucleation processes and near 
bifurcation points they can play an essential role. 

More importantly, the above calculations and estimates show that under not very extreme condi- 
tions, even less extreme if diffusion coefficients are smaller than our typical value, one may observe 
the consequences of the interactions between internal fluctuations and periodic or chaotic dynamics 
in the spatial domain. 
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VI. CONCLUSION 



The reactive lattice-gas automaton provides a description of the system at a mesoscopic level. 
On the smallest scales one has individual reactive events and random-walk dynamics. Of course, 
in general, at this level one must regard the particles as fictitious, representing some collection 
of real particles. However, one may consider physical circumstances where the particle numbers 
in the automaton correspond closely to those of a real system. In this case, the dynamics may 
be considered to be a model of real collision dynamics. At the coarse-grained level one observes 
a noisy version of the periodic or chaotic attractor. At this mesoscopic level one has interplay 
between the local periodic or chaotic dynamics, internal fluctuations and diffusion to give rise to 
the observed spatiotemporal structure. 

In the chaotic regime internal fluctuations are amplified by the dynamcis and give rise to am- 
plitude fluctuations in the coarse-grained concentration variables. Provided the diffusion length is 
large enough the system maintains phase coherence over all space and the phase space trajectory 
of the globalUy-averaged concentrations lies on a noisy periodic attractor. In the periodic regime 
fluctuations from the deterministic limit cycle exhibit linear growth and a noisy limit cycle is also 
observed in this periodic regime. This behavior is consistent with earlier considerations of the 
noisy dynamcis of coupled map lattices. jl9j If the diffusion length is sufficiently small then the 
chaotic system will no longer be able to maintain phase coherence and phase turbulence will ensue, 
similar to that seen in oscillatory media. 

The reactive lattice-gas automaton is a synchronously-updated, probabilistic cellular automaton 
and thus shares some of the features that were seen to be essential in the observation of non-trivial 
global attractors. However, the well-stirred dynamics for large systems was seen to closely 

approximate that of the continuous-time mass action ODE. One may show in this limit that the 
mean field automaton equations are simply an Euler discretization of the mass action rate law. 
The time step in this discretization is controlled by the scale factor h introduced in ( |ll] ) and its 
small size insured that the automaton dynamics in this limit yielded results comparable to the 
solution of the underlying ODE. The reative lattice-gas model is intrinsically statistical in nature 
and encorporates the effects of internal fiuctuations. Thus, there are a number of subtle issues 
involving the passage between disctete-time and continuous-time dynamics, and their implications 
for systems with broken time translation symmetry, that merit further investigation in the context 
of these models. 
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FIG. 1. Bifurcation diagram constructed from the solution of the mass-action equations (2). Here u and 
w, respectively, are 17 and W concentrations on a Poincare section P, defined as {P : Pv = Ve, ^pu, Pw > 0}. 
The bifurcation diagram shows the projection onto u in the Poincare section. 
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FIG. 3. 3-D phase space plots of the global concentrations for K2 = 11.52 for various system sizes: (a) 
64 X 64 lattice, (b) 128 x 128 lattice, (c) 256 x 256 lattice and (d) 512 x 512 lattice. 




FIG. 4. The twelve panels of this figure show the local concentrations of species U for various times 
indicated in Fig. |^. Time increases from left to right and top to bottom and the total time spanned is 
about one and a half cycles on the attractor. The concentration is coded by gray shades with black high 
concentration and white low concentration of U. 
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FIG. 5. Plot of the global concentration of species U in the time interval corresponding to the panels 
showing the spatial structure of the the U concentration field. The points on the curve are taken at equally 
spaced time intervals. The panels in Fig. ^ correspond to the open squares while the open triangles denote 
intermediate times where the structure is not displayed. 



FIG. 6. The six panels of this figure show the local concentrations of species f/ on a 512 x 512 lattice for 
various times within one cycle on the global attractor. Time increases from left to right and top to bottom. 
The concentration is coded by gray shades with black high concentration and white low concentration of 
U. 
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FIG. 7. Plot of the coaxse-grained concentration of species U, Pu(r, t), in four cells as a function of time 
for K2 = 11.52. 
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FIG. 8. Global phase space dynamics for K2 = 12.8, a parameter value in the periodic regime, for (a) 
the deterministic mass action rate law (B) and (b) automaton dynamics on a 128 x 128 lattice. 
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FIG. 9. Plot of the coarse-grained concentration of species U, p„(f, i), in four cells as a function of time 
for K2 = 12.8. 
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FIG. 10. Phase space plots ol the global atractors determined from automaton dynamics on 128 x 128 
lattices for two values of the diffusion coefficient: (a) Dt = D/6 and (b) Dt = D/10. 
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FIG. 11. Plot of the coarse-grained concentration of species U, pu{r,t), in four cells as a function of 
time for K2 = 11.52 for = D/6. 
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